On the structure of correlations in the three dimensional spin glasses. 
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We investigate the low temperature phase of three-dimensional Edwards- Anderson model with 
Bernoulli random couplings. We show that at a fixed value Q of the overlap the model fulfills the 
clustering property: the connected correlation functions between two local overlaps decay as a power 
whose exponent is independent of Q for all < \Q\ < Qea- Our findings are in agreement with the 
RSB theory and show that the overlap is a good order parameter. 



Spin glasses have unusual statistical properties. In 
mean field theory there are intensive quantities that fluc- 
tuate also in the thermodynamic limit. This is the effect 
of the coexistence of many equilibrium states. The corre- 
lations functions inside a given state should have a power 
law behavior: below the critical temperature spin glasses 
are always in a critical state (many glassy systems should 
share this behavior) . These predictions of mean field the- 
ory have been never studied in details, apart from [tJ; the 
aim of this paper is to address this point in a careful way. 

In order to better characterize the behavior of spin 
glasses, it is convenient to consider two replicas, or clones, 
of the same system (let us call them a{i) and r(i) where 
i denotes the point of the lattice). The two clones share 
the same Hamiltonian Hj, the label J indicates the set 
of random coupling constants. 

Let us define the local overlap q{i) = (j{i)T{i) and the 
global overlap q = J^iQi'^)^ ^ being thejotal vol- 
ume. For the three dimensional EA model 3] at zero 
magnetic field (to be defined later) all simulations con- 
firm that the probability distribution of Pj (q) is non triv- 
ial in the thermodynamic limit, it changes from system to 
system and its average over the disorder, that we denote 
as P{q) = E[Pj(q)], is non trivial and it has a support in 
the region from —qsA to qsA, qsA being the overlap of 
two generic configurations belonging to the same state. It 
is usually assumed that the function P{q) has in the infi- 
nite volume limit a delta function singularity at q = qea, 
that appears as a peak in finite volume systems. In the 
presence of multiple states the most straightforward ap- 
proach consists in identifying the clustering states (i.e. 
those where the connected correlation functions go to 
zero at large distance) and to introduce an order para- 
meter that identifies the different states. This task is 
extremely difficult in a random system where the struc- 
ture of the states depends on the instance of the system. 
However the replica theory is able to make predictions 
without finding out explicitly the set of states. At this 
end the introduction of the two clones plays a crucial role. 
Indeed if the global overlap q has a preassigned value, the 
theory predicts that the correlations of local overlaps q{i) 
go to zero at large distances. In other words g is a good 



order parameter. 

For each realization of the system we consider two 
clones. The observables are the local overlaps q{i) and 
their correlations. We define {0)q the expectation value 
of the observable O in the J-dependent Gibbs ensemble 
restricted to those configurations of the two clones that 
have global overlap q = Q. We define the average expec- 
tations values {0)q as weighted average over the systems 
of restricted expectation values [3| : 



{0)q = 



E[P.jiQ){OyQ] 

E[pjm 



(1) 



The main statement of the RSB theory [l| (to which 
we bring evidence in this paper) is that the Q-dependent 
connected correlation functions go to zero when com- 
puted in the ensemble (■)q, i.e. the states (•)q are clus- 
tering. The procedure is very similar to the one used in 
ferromagnetic models to construct clustering states (for 
example in the Ising case by considering averages with 
positive, or negative, total magnetization). The overlap 
constraint state is not clustering when the equilibrium 
state is locally unique (apart from a global change of 
signs) Q. Indeed this is the only known example of a 
system where the clustering states are labeled by a con- 
tinuously changing order parameter in absence of a con- 
tinuous symmetry, like rotations or translations. 

This clustering property has far reaching consequences: 
for example the probability distribution P{q)^ of the 
window overlaps i.e. the average overlap over a re- 
gion of size W, becomes a delta function in the infinite 
volume limit: limn/_,oo P{q)^ = ^{Q ~ o)- Indeed when 
the region W becomes large, the window overlaps are in- 
tensive quantities that do not fluctuate inside clustering 
states. 

It is convenient to recall the known theoretical results 
for the connected correlation functions in the case of 
short range Ising spin glasses: 



Gix\Q) = {qix)q{0))c 



C{x\Q) ^ G{x\Q) - (2) 



and their Fourier transforms C{k\Q). 
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Figure 1: The correlation function at distance 10 G{10\Q) 
averaged in 20 bins (round points) and in ten bins (crosses) 
of Q^, versus the average of inside the bin. The data are 
for a system of size 20 and the straight is the best fit to the 
data. 



Figure 2: Circles are the value of the exponent a{Q) by fitting 
x'°'(Q)i as a power of L: for each value of we show two 
points corresponding to s = 1 and to s = 2. Squares repre- 
sent the same quantity a{Q) as obtained through the scaling 
approach visible in fig. |3l 



The simplest predictions are obtained starting from 
mean field theory and computing the first non trivial 
term Q Neglecting logarithms we have in the small k 
region 



C{k\Q) oc k-^ for 

C{k\Q) oc for 

C{k\Q) cx k-^ for 

C(fc|Q)oc(fc2+^(g)-2)-i for 



g = o , 

< Q < qea , 

Q = QEA , 

Q > qsA ■ (3) 



These results are supposed to be exact at large distances 
in sufficiently high dimensions, i.e. for D > 6. 

The reader may be surprised to find a result for Q > 
Qea because the function P{q) is zero in this region in 
the infinite volume limit. However for finite systems P{q) 
is different from zero for any q albeit it is very small 
lol . 11 1 in the region q » qea- For Q > qEA an analytic 
computation of the function C{k\Q) has not yet been 
done, however it is reasonable that the leading singularity 
near to fc = in the complex plane is a single pole, leading 
to an exponentially decaying correlation function [ISj] . 

When the dimensions become smaller than 6 we can 
rely on the perturbative expansion in e = 6 — D, where 
only the first order is (partially) known [T^j. It seems 
that predictions at Q = qea should not change and the 
form of the k = singularity at Q = qEA (i-e. when 
the two clones belong to the same state) remains k^^ 
as for Goldstone Bosons. On the contrary the fc = 
singularities at Q < qEA should change. We can expect 
that 



C{k\Q) cx for < Q < 



qEA 



(4) 



These perturbative results are the only information we 
have on the form of a{Q). The simplest scenario would 



be that a{Q) is discontinuous at Q = and constant in 
the region < Q < qEA- It would be fair to say that 
there is no strong theoretical evidence for the constancy 
of a{Q) in the region < Q < qsA, apart from generic 
universality arguments. On the contrary the discontinu- 
ity at Q = could persist in dimensions not too smaller 
than 6, and disappear at lower dimensions, as supported 
by our data in Z? = 3. 

In three dimensional case in configuration space we 
should have: 



C{x\Q) oc x""^'^^ for < g < 



qEA 



(5) 



with a{qEA) = 1 (indeed in general we have that a{Q) = 
D — a{Q)). For Q > qEA the correlation should go to 
zero faster than a power: we tentatively assume that 



C(x\Q)(xx ^ exp(-x/^((5)) for qEA < 



(6) 



(also other behaviours are possible). 

In this paper we will numerically study the properties 
of the two overlaps connected correlations functions in 
the three dimensional EA model. The Hamiltonian of the 
EA model [3|] is given by — — jfiajwith 
Jij = ±1 (symmetrically distributed) and Ising spins 

(T./-±l. 

We have studied cubic lattice systems with pe- 
riodic boundary conditions of side L for L = 
4, 6, 8, 10, 12, 16, 20. The simulation parameters are the 
same used in reference Q. We present the results only 
at temperature T = 0.7, while the critical temperature is 
about T = 1.11. 

We have first classified the configurations created dur- 
ing the numerical simulations according to the value of 
the global overlap q. Since the properties of the configu- 
rations are invariant under the symmetry (q ~<z), we 
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Figure 3: The quantity L°'^'^^ CL{x\Q)g{x / L) with g{z) = 
{l/z + 1/(1 — versus sin(7rx/L) for the set of data in 

the first bin < 0.5, using the best value of a(Q). The left 
panel display the data for all the correlation at distances x>l 
(the corresponding value of a{Q) being 0.50); in the right 
panel we have only the data with x > 2 (the corresponding 
value of a(Q) being 0.54). 



have classified the configurations into 20 equidistant bins 
in q^: e.g. the first bin contains all the configurations 
where < < 1/20. In this way we compute the corre- 
lations C{x\Q). We have measured the correlations only 
along the axes of the lattice: x is an integer restricted 
to the range 0, L/2. As a control we have done the same 
operation with 10 bins obtaining similar results. 

We have firstly verified that the connected correlations 
vanish for large systems. At this end in fig. ([T]) we have 
plotted for L — 20 (our largest system) the correlation 
G{10\Q) versus the average of in the bin. We see 
that the two quantities coincide. The data show a strong 
evidence for the vanishing of the connected two point 
correlation function. The predictions of replica theory is 
G'(10|Q) = Q^, neglecting small corrections going to zero 
with the volume. 

Further information can be extracted from the data. 
The analysis of the data should be done in a different 
way in the two regions < < q'^^ and q^^ < as 
far as two different behaviour are expected. In our case 
can be estimated to be around 0.4. 

In the region < < q'^^^ the power law decrease 
([5]) of the correlation is expected. To test this hypothesis 
[3] we define for each L the quantities x'l\Q) 



L/2 



(7) 



where Cl{x\Q) is the connected correlation function in 
a system of size L, i.e. Gl{x\Q) — [lB|- For large L 
Xl\Q) should behave as L''+'^~°'(Q\ We have evaluated 
the previous quantity [17l | for s = 1, 2. In the region 



Q2 < 0.4 we have found that the ratio Xl\Q)/Xl\Q) 
is well linear in L. Here the data for Xl ^ (Q) can be well 
fitted as a power of L and the exponents a{Q) computed 
using s = 1 and s = 2 coincide within their errors. These 
results are no more true in the region 0.5 < indicating 
that a power law decrease of the correlation is not valid 
there. The exponents we find with this method are shown 
in fig. dH). 

In order to check these results for a{Q) we have used 
a different approach. In the large volume limit the cor- 
relation function should satisfy the scaling 



(8) 



The value of a{Q) can be found by imposing this scal- 
ing. At this end for each value of Q we have plotted 
ia(Q)(7^(a;|Q) and found the value of a{Q) for which 
we get the best collapse. The result of the collapse is 
shown in fig. ([3]) for Q around zero, where for graph- 
ical purpose we have plotted L"'-Q^CL{x\Q)g{x/L) ver- 
sus sin(7rx/L), where the function g has been added to 
compress the vertical scale (we find convenient to use 
g{z) = (1/z -I- 1/(1 - z))-"^?), following 0). In the left 
panel we show the collapse using all points with a; > 1, 
in the right panel we exclude the correlations at distance 
X — 1. The corresponding values of the exponent are 
shown in fig. ^ and they agree with the ones coming 
from the previous analysis in the region of < 0.4. 

The exponent a{Q) is a smooth function of which 
goes to 1 near — 0.4 in very good agreement with 
the theoretical expectations. We do not see any sign 
of a discontinuity at Q = 0, and this is confirmed by an 
analysis with an high number of bins (e.g. 100). However 
it is clear that for lattice of this size value we cannot 
expect to have a very high resolution on Q and we should 
look to much larger lattices in order to see if there is a sign 
of a building up of a discontinuity and of a plateaux. The 
value of the exponent that we find at Q = is consistent 
with the value 0.4 found from the dynamics [l^, and 
with the value 0.4 found with the analysis of the ground 
states with diflFerent boundary conditions 0. 

From the previous analysis it is not clear if the expo- 
nent a{Q) has a weak dependence on Q or if the weak 
dependence on Q is just a pre- asymptotic effect. In or- 
der to clarify the situation it is better to look to the 
connected correlations themselves. In fig. (jlj we display 
the connected correlation Cl{x\Q) as function of for 
X = 6, 7, 8, 9 at L = 20, our largest lattice (for the result 
at X = 1 see [HI). We can fit the correlations at fixed L 
(e.g. L = 20) for large x as: 

Cl{x\Q)^A{x,L){Q^-B{x,L)^), Q^<B{x,L)^ (9) 

while Cl{x\Q) is very near to zero for > B{x,L)^. 
The goodness of these fits improves with the distance 
(similar results are valid at smaller L). The value of 
B{x, LY is near to q^j,^ and it is slightly decreasing with 
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Figure 4: The connected correlation Cl{x\Q) as function of 
for x = 6, 7, 8, 9 at L = 20. The straight lines are linear 

fits. 



Figure 5: Correlations functions at L = 20 for — 0.475, 
0.525, 0.575, 0.625, versus y = (1/x + 1/(2L - x)y\ The 
straight line is proportional to ■ 



L. The validity of the fits Q for large L would imply 
that in the region \q\ < qea the large distance decrease 
of the correlation function should be of the form A{Q^ — 
q%ji}f{x) and therefore the exponent a{Q) should not 
depend on Q. 

However near q = qea we should have a real crossover 
region. In Fig. \5\we show CL{x\Q)g{x) at L = 20 for 
0.475 < Q2 < 0.625, versus y = {1/x + 1/(2L - x))-^ 
(we use the variable y = x{l — 0{x/L)) to take care 
of finite volume effects) with 17(2;) = (1 — 2x/L)~'^. It 
seems that the data at > 0.475 decrease faster than 
a power at large distances and that the data at = 
0.475 are compatible with a power with exponent — 1 
It is difficult to extract precise quantitative conclusions, 
without a careful analysis of the L dependence. We hope 
that this will be done when the data on the correlation 
functions on larger lattices will be available. 

In the region < our task is different: the cor- 
relations are short range and we would like to compute 
if possible the correlation length. At this end we have 
fitted the data as 

Cl{x\Q) — exp(— a;/^i((5)) + (x ^ L — x) + const 

X + 1 

(10) 

The choice of the fit is somewhat arbitrary, however we 
use it only to check that the correlation length diverges 
at qEA and that near qEA are well fitted by a l/x power. 
The fits are good, but this may not imply the correctness 
of the functional form in eq. (fTU)) . We find that far from 
Q = 0.5 the correlation length is independent of L, (it 
is quite small). We have tried to collapse the data for 
L > 8 in the form 

UQ)^Lf{{Q'~ql^)L^/n ■ (11) 
A reasonable collapse has been obtained, however the 



qEA is quite small (i.e. 0.25): it is quite possible that 
there are finite volume effects and different ways to eval- 
uate qEA give different results on a finite lattice and they 
would converge to the same value in the infinite volume 
limit. 

In conclusion the global overlap for a two clone system 
is a well defined order parameter such that in the appro- 
priate restricted ensemble the two points connected cor- 
relation function decays at large distance. The connected 
correlations decay as a power whose exponent seem to be 
independent from Q for < \Q\ < qsA- the value of the 
exponent is in agreement with the results obtained in 
different context at Q = 0. Moreover the connected two 
points correlation functions ai Q = qEA decays like 1/x 
in agreement with the detailed predictions coming from 
replica theory. 
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